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(57) Abstract: The invention concerns a method for fine synchronisation 
on a reception signal (S) corresponding to a reference signal (TS) trans- 
mitted in a transmission channel, comprising the following steps: select- 
ing a signal source producing a characterisation signal (X) following its 
passage in the transmission channel; establishing a characterisation ma- 
trix (L) to estimate the characterisation signal (X) co variance; identifying 
characteristic modal values which are the highest characteristic values (Xj) 
of said characterisation matrix (L); computing the source signal correla- 
tion function c(t) with the sum of characteristic vectors (v$ associated 
with the characteristic modal values; searching said correlation function 
c(t) first maximum. 



OS 

O 



Xi i c i 



9 ' • • f 



Vi 1 a 1, d 



g(t) 





d 




z g<t). v. 




i = l 



Max[c(t)] -> t 0 



(57) Abrege: L* invention concerne un procede de synchronisation fine 
sur un signal de reception (S) correspondant a un signal de reference (TS) 
emis dans un canal de transmission. II comprend les Ctapes suivantes: se- 
lection d*un signal source produisant un signal de caracterisation (X) suite 
a son passage dans le canal de transmission, etablissement d'une matrice 
de caracterisation (L) pour estimer la covariance du signal de caracterisa- 
tion (X), identification des valours propres dominances qui sont les valeurs 
propres (X*) les plus elevees de cette matrice de caracterisation (L), calcul 
de la fonction de correlation c(t) du signal source avec la somme des vec- 
teurs propres (Vj) associes aux valeurs propres dominantes, recherche du 
premier maximum de cette fonction de correlation c(t). 
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The invention concerns a method for fine synchronisation on a 
reception signal (S) corresponding to a reference signal (TS) transmitted in 
a transmission channel, comprising the following steps: selecting a signal 
source producing a characterisation signal (X) following its passage in the 
transmission channel; establishing a characterisation matrix (L) to estimate 
the characterisation signal (X) covariance; identifying characteristic modal 
values which are the highest characteristic values (AO of said characterisation 
matrix (L); computing the source signal correlation function c(t) with the sum 
of characteristic vectors (vj) associated with the characteristic modal values; 
searching said correlation function c(t) first maximum. 

(57) Abrege 

L* invention concerne un proc6d6 de synchronisation fine sur un signal 
de reception (S) correspondant a un signal de r6f6rence (TS) ernis dans un 
canal de transmission. II comprend les 6tapes suivantes: selection d'un signal 
source produisant un signal de caract6risation (X) suite a son passage dans 
le canal de transmission, dtablissement d'une matrice de caract6risation (L) 
pourestimer la covariance du signal de caracterisation (X), identification des 
valeurs propres dominantes qui sont les valeurs propres (A;) les plus elevees 
de cette matrice de caractdrisation (L), calcul de la fonction de corr61ation 
c(t) du signal source avec la somme des vecteurs propres (vj) associds aux 
valeurs propres dominantes, recherche du premier maximum de cette fonction 
de correlation c(t). 
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Proced~ de synchronisation fine sur un signal regu 
d'un canal de transmission 
La presente invention concerne un procede de 
synchronisation fine sur un signal de reception 
5 correspondant a un signal de reference emis dans un canal de 
transmission . 

Dans un systeme de transmission, notamment par ondes 
radios, un emetteur emet un signal de reference dans un 
canal de transmission a destination d'un recepteur. Une des 

10 premieres operations que doit realiser le recepteur est la 
synchronisation sur le signal de reception. Ce probleme est 
bien connu de l'homme du metier et il n'est done pas 
necessaire de rappeler ici les differentes techniques qui 
sont employees pour obtenir cette synchronisation. 

15 Lorsque les signaux sont numeriques, on a coutume 

d'evaluer l'ecart de synchronisation au moyen d'une unite de 
temps, le temps bit, qui est l'ecart temporel separant deux 
bits successifs d'un signal. Or, il apparait que les 
solutions disponibles dans l'etat de l'art ne permettent pas 

2 0 d'acquerir une synchronisation beaucoup plus precise que le 
temps bit. 

Cette precision peut s'averer insuffisante dans 
certains cas. En effet, e'est la synchronisation qui donne 
le temps de trajet du signal dans le canal de transmission 

25 ou, autrement dit, le temps de transmission entre 1' emetteur 
et le recepteur. Cette donnee est importante car, dans un 
systeme de radiocommunication duplex oCi une station de base 
est en communication avec un terminal , chaque equipement 
etant pourvu d'un emetteur-recepteur , le terminal doit 

30 fonctionner de sorte que le signal qu'il emet arrive a un 
instant precis en reference a la base de temps de la station 
de base. Pour ce faire, il faut naturellement que le 
terminal connaisse le temps que met ce signal pour parvenir 
a la station de base. 

35 D 1 autre part, ce temps de transmission reflate 

directement la distance separant l 1 emetteur du recepteur. On 
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comprend bien que la precision de cette distance est 
f ondamentale lorsqu'il s'agit de proceder a la localisation 
du terminal en reperant sa position relative par rapport a 
une ou plusieurs stations de base. Le probleme general de la 
localisation d'un terminal est une preoccupation tres 
actuelle du fait de ses applications au nombre desquelles on 
citera, par exemple, les strategies de changements de 
cellules ("handover") dans les reseaux cellulaires. On 
mentionnera egalement le domaine de la s6curite, qu'il 
faille situer geographiquement la provenance d'un appel 
d'urgence ou bien la position d'un vehicule vole equipe du 
terminal . 

La presente invention a ainsi pour objet un procede de 
synchronisation dont la precision est bien superieure au 
temps bit. 

Selon 1' invention, le procede de synchronisation fine 
sur un signal de reception correspondant a un signal de 
reference emis dans un canal de transmission comprend les 
6tapes suivantes : 

- selection d'un signal source produisant un signal de 
caracterisation suite a son passage dans le canal de 
transmission , 

- etablissement d'une matrice de caracterisation pour 
estimer la covariance du signal de caracterisation, 

- identification des valeurs propres dominantes qui 
sont les valeurs propres les plus elevees de cette matrice 
de caracterisation, 

- calcul de la fonction de correlation du signal 
source avec la sorame des vecteurs propres associes aux 
valeurs propres dominantes, 

- recherche du premier maximum de cette fonction de 
correlation. 

Lorsque l 1 increment de temps adopte pour le calcul de 
la fonction de correlation est choisi suffisamment petit, en 
tout cas bien inferieur a un temps bit, ce procede permet 
d'obtenir une tres bonne precision. 
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Selon une premiere option, le nombre de ces valeurs 
propres dominantes est predetermine. Typiquement, ce nombre 
represente de l'ordre de 20 a 30 % de la dimension de la 
matrice de caracterisation. 
5 Selon une deuxieme option , le rapport de la somme des 

valeurs propres dominantes a la somme de toutes les valeurs 
propres est superieur ou egal a un nombre predetermine, Ici, 
le nombre retenu sera souvent superieur a 90 %, 95 % par 
exemple. 

10 Selon une troisieme option, le procede comprenant de 

plus une etape d' estimation du bruit additif dans le canal 
de transmission, les valeurs propres dominantes sont telles 
que leur somme soit inferieure ou egale a la somme de toutes 
les valeurs propres diminuee du bruit additif. 

15 De plus, l 1 estimation du bruit additif est realisee en 

normalisant le bruit instantane qui est evalue au moyen du 
signal de reception, du signal de reference et d'une 
estimation de la reponse impulsionnelle du canal de 
transmission . 

20 Avantageusement , en notant A la matrice de 

transmission associee au signal de reference, l 1 expression 
du bruit instantane est la suivante : Nq = S - A.X. 

Par ailleurs, quelle que soit l 1 option eventuellement 
retenue, la matrice de caracterisation resulte d'une 
25 operation de lissage. 

Selon un mode de realisation pref erentiel , le signal 
de caracterisation est une estimation de la reponse 
impulsionnelle du canal de transmission. 

On peut egalement prevoir que le signal de 
30 caracterisation soit le signal de reception. 

La presente invention apparaitra maintenant de maniere 
plus detaillee dans le cadre de la description qui suit ou 
sont proposes des exemples de mise en oeuvre a titre 
illustratif, ceci en reference aux figures annexees qui 
35 representent : 
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- la figure 1, une premiere variante de realisation de 
1 1 invention , et 

- la figure 2, une deuxieme variante. ' A 

Les elements communs aux deux figures sont affectees 
5 d'une seule et meme reference. 

Le recepteur a deja acquis une synchronisation 
grossiere sur le signal de reception, de 1'ordre du temps 
bit, au moyen de l'une quelconque des solutions disponibles* 
Ce signal de reception correspond a un signal de 
10 reference produit par l'emetteur et connu du recepteur* Ce 
signal de reference peut etre connu a priori, c'est-a-dire 
qu'il s'agit d'une sequence d 'apprentissage formee de 
symboles identifies. II peut egalement etre connu a 
posteriori au moyen de techniques generiquement referencees 
15 sous le terme de sondage aveugle. Dans ce cas, au cours de 
la procedure de synchronisation, le recepteur regenere la 
suite des symboles formant le signal de reference a partir 
du signal de reception. 

II convient en premier lieu de caracteriser la 
20 transmission entre l'emetteur et le recepteur* A cet effet, 
on selectionne un signal source qui, produit par l'emetteur, 
donne au niveau du recepteur, apres transmission dans le 
canal, un signal de caracterisation. 

Naturellement, si le signal source est le signal de 
25 reference, ce signal de caracterisation est le signal de 
reception lui-meme. Ce n f est pas cependant pas tou jours la 
solution optimale quant a la complexite et aux performances 
du procede de 1 1 invention . 

Une autre solution consiste a retenir une impulsion 
30 modulee comme signal source, le signal de caracterisation 
devenant alors la reponse impulsionnelle du canal de 
transmission . > 

A titre d'exemple, le systeme de radiocommunication 
cellulaire numerique GSM fait appel a une sequence 
35 d'apprentissage de 26 symboles, la reponse impulsionnelle 
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etant generalement estimee avec 5 coefficients puisqu'on 
admet que la dispersion du canal vaut 4 . 

Dans ce cas, le signal de reception a une dimension 
maximale de 22 qui est sensiblement plus importante que 
5 celle de la reponse impulsionnelle. 

On examinera done successivement deux variantes de 
mise en oeuvre de 1 ■ invention en commengant par le cas ou le 
signal source est une impulsion modulee en reference a la 
figure 1 . 

10 L 1 estimation de la reponse impulsionnelle ne pose pas 

de difficultes en soi car de nombreuses methodes permettent 
d'y parvenir, par exemple la methode dite du critere des 
moindres carres qui est decrite notamment dans les demandes 
de brevet FR 2696604 et EP 0564849. En matiere de rappel , 

15 cette technique fait appel a une matrice de mesure A 
construite a partir de la sequence d 1 apprentissage TS de 
longueur n. Cette matrice comprend (n-d) lignes et (d+1) 
colonnes, d representant la dispersion du canal. L 1 element 
figurant a la ieme ligne et a la jeme colonne est le (d+i- 

20 j)ieme symbole de la sequence d 1 apprentissage, soit en 
notant a^ le ieme symbole d'une sequence TS de 26 symboles : 

r a4 a3 a2 ai ao ^ 
as a4 a3 a2 ai 
a6 as a4 a3 a2 



\^25 &21 J 

La sequence d 1 apprentissage est choisie telle que la 
25 matrice A^A soit inversible ou l'operateur represente la 
transposition . 

Dans le signal de reception S, on ne prend pas en 
compte les quatre premiers symboles sq a S3 car ceux-ci 
dependent egalement de symbole inconnus emis avant la 
30 sequence d 1 apprentissage, etant donne que la dispersion du 
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canal vaut 4. Par un abus de langage on definira done 
dorenavant le signal de reception comme un vecteur S ayant 
pour composantes les symboles regus, S4, S5, . .. , S25. 

Des lors, 1 1 estimation de la reponse impulsionnelle X 
5 prend la forme suivante : 

X = (A fc A)" 1 A fc . S 
L'etape suivante du procede de 1' invention consiste a 
etablir une statistique de cette reponse impulsionnelle. Par 
statistique, on entend un ensemble de donn§es refletant la 
10 valeur moyenne de cette reponse sur une p6riode d 1 analyse. 

On construit done une matrice de lissage des 
differentes estimations X obtenues pendant la periode 
d' analyse pour obtenir une estimation de la covariance 
associee a cette reponse impulsionnelle. On entend ici 
15 lissage dans un sens tres general, e'est-a-dire toute 
operation permettant de lisser ou de moyenner la reponse 
impulsionnelle sur la periode d' analyse. 

Un premier exemple de lissage consiste a effectuer la 
moyenne de la matrice XX h sur la periode d' analyse supposee 
20 comprendre m sequences d'apprentissage : 

L(XX h ) = ^-^ XX h 

L'operateur .h represente la transformation 
hermitienne ou transposition conjugaison complexe. 
25 Un second exemple de lissage consiste a actualiser, 

apres reception de la ieme sequence d'apprentissage, la 
matrice de lissage Li-^XX* 1 ) obtenue a la (i-l)ieme sequence 
d'apprentissage au moyen d'un coefficient multiplicatif a, 
ce facteur etant generalement connu sous le nom de facteur 
30 d'oubli de lissage et etant compris entre 0 et 1 : 

Li(XX h ) = aXiXi h + (l-cOLj-ifXxh) 
L 1 initialisation peut se faire par tous moyens, 
notamment au moyen de la premiere estimation X obtenue ou 
bien par une moyenne obtenue comme ci-dessus pour un faible 
35 nombre de sequences d'apprentissage. 
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Dans un souci de simplification, la matrice de lissage 
L(XX tl ) qui est en fait une matrice de caracterisation 
statistique, sera desormais notee L. 

Le procede comprend ensuite une etape de recherche des 
5 couples (valeur propre, vecteur propre) de la matrice de 
caracterisation . 

Cette etape ne sera pas plus detaillee car bien connue 
de l'homme de metier. 

Les valeurs propres X± sont maintenant classees, par 
10 ordre decroissant. En effet, la somme de ces valeurs 
correspond a l 1 energie du signal de caracterisation X 
composee pour partie d'un signal utile qui est I 1 image du 
signal source et pour partie du bruit additif N du canal de 
transmission * 

15 II vient que les valeurs propres dominantes, celles 

qui sont le plus elevees, representent le signal utile, 
tandis que les valeurs propres les plus faibles representent 
le bruit. 

Selon une premiere option, le procede consiste a 
20 retenir un nombre predetermine de valeurs propres 

dominantes. Par exemple, pour une reponse impulsionnelle a 

5 coefficients, on retient les deux premieres valeurs 

propres X± et X2* 

Selon une deuxieme option, on considere que le signal 
25 utile presente une energie qui est une fraction 

predeterminee f de l 1 energie du signal de caracterisation. 

Ainsi en notant X± les valeurs propres pour i variant de 1 a 

p, il y aura d valeurs propres dominantes, d etant obtenu 

comme suit : 



30 



d d + 1 

P P 

1=1 1 = 1 

La fraction f peut §tre fixee a priori, a une valeur 
de 95 % par exemple. Cette fraction peut egalement deriver 
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du rapport signal a bruit du signal de reception obtenu par 
ailleurs. 

Selon une troisieme option du procede, sans doute la 
plus performante, le bruit additif N est estime directement 
5 a partir du signal de reception et de la matrice de mesure 
A. En effet, en notant Nq le vecteur bruit affectant le 
signal de reception, il vient que : 
S = AX + N 0 

Compte tenu du fait que les vecteur s S et Nq ont 22 
10 composantes, le bruit additif N peut s'exprimer de la 
maniere suivante : 

N = (^) (S - AX) h (S - AX) 

Naturellement, cette estimation du bruit additif peut 
etre moyennee en lissee. 
15 En reprenant les notations precedentes et en 

normalisant les energies, il vient que : 

P P 

E X. > N et £ A,* < N 

i=d i=d+l 
On obtient done les valeurs propres dominantes a 

20 partir d'une estimation directe du bruit* 

Quelle que soit l 1 option precedemment retenue, l f etape 
suivante du procede consiste a calculer la fonction de 
correlation du signal source avec la somme des vecteurs 
propres v^ associes aux valeurs propres dominantes X^. 

25 Le signal source est surechantillone par rapport au 

temps bit et on le notera done g(t) ou t qui represente le 
temps est une variable discrete dont le pas de 
quantification vaut, a titre d'exemple, 1/32 temps bit. Il 
est represente par un vecteur de meme dimension que le 

30 signal de caracterisation, soit 5 dans I'exemple adopte. La 
fonction de correlation c(t) est calculee par exemple pour t 
variant de -1 6 +1 temps bit au moyen de !• expression 
suivante : 
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c(t) = z g(t). v. 
i = l 

Le point figurant entre le signal source g(t) et le 
vecteur propre v^ represente classiquement le produit 
scalaire. 

5 La derniere etape du procede consiste a rechercher la 

valeur to de t la plus proche de zero qui correspond au 
premier maximum relatif de la fonction correlation c(t). 
C'est cette valeur particuliere t 0 qui donne I'ecart de 
synchronisation recherchee par rapport au signal de 
10 reception. 

Par ailleurs, on peut considerer la fonction 
complementaire c'(t) suivante : 

c f (t)= z g(t).v 

i = d + 1 

II faut remarquer que la valeur particuliere to 
15 mentionnee ci-dessus peut egalement etre obtenue en 
recherchant la valeur de t la plus proche de zero qui 
correspond au premier minimum relatif de la formation 
complementaire c ■ ( t ) . 

Ces deux methodes pour obtenir I'ecart de 
20 synchronisation t 0 sont done equivalentes . 

En refeence a la figure 2, considerons maintenant une 
deuxieme variante de 1' invention selon laquelle le signal de 
caracterisation est le signal de reception S, si bien que le 
signal source est maintenant le signal de reference, soit 
25 dans le cas du GSM, la sequence d 'apprentissage TS modulee 
GMSK (pour "Gaussian Minimum Shift Keying"). 

La statistique du signal de caracterisation est done 
estimee au moyen d'une matrice de caracterisation qui est 
maintenant obtenue par lissage des differentes occurrences 
30 du signal de reception S. Ici encore, on considere le terme 
lissage dans un sens tres general. 

La matrice de caracterisation L prend done la forme 
suivante : 
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H 1 111 K 
L(SS n ) = — Y SS n ou bien, 
m 

Li(SSh) = aSiSih + (l-cOLi^SS 11 ) 

On recherche ensuite les p 1 valeurs propres A,' ^ de 
cette matrice et, comme dans la premiere variante, on 
identifie les d' valeurs propres dominantes. 

On calcule maintenant la fonction de correlation f(t) 
de la sequence d 1 apprentissage modulee et de la somme des 
vecteurs propres v'^ associes aux valeurs propres dominantes 
^* ' i * 

La encore, la sequence d ' apprentissage g'(t) est 
surechantillonnee et elle est representee par un vecteur de 
2 2 composantes. La fonction de correlation devient done : 

d 1 

f(t) = z g f (t).v. 

i = 1 1 

Comme auparavant, on peut definir une nouvelle 
fonction complementaire f'(t) : 

P' 

f f (t)= Z g'(t).v'. 
i = d ■ + l 

Le proc§de se termine de la meme maniere en 
recherchant le premier maximum de la fonction de correlation 
f(t) ou le premier minimum de la fonction complementaire 
f '<t). 

L 1 invention peut ainsi etre mise en oeuvre de 
differentes manieres, le point essentiel etant de disposer 
d'un signal source et du resultat de sa transmission, a 
savoir du signal de caracterisation . 
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RE VEND I CAT I ONS 

l)Procede de synchronisation fine sur un signal de reception 
(S) correspondant a un signal de reference (TS) emis dans un 
canal de transmission, caracterise en ce qu'il comprend les 
5 etapes suivantes : 

- selection d'un signal source produisant un signal de 
caracterisation (X, S) suite a son passage dans ledit canal 
de transmission, 

- etablissement d'une matrice de caracterisation (L) pour 
10 estimer la covariance dudit signal de caracterisation (X, 

S), 

- identification des valeurs propres dominantes qui sont les 
valeurs propres (X±, X % jj les plus elevees de cette matrice 
de caracterisation (L), 

15 - calcul de la f one t ion de correlation (c(t), f(t)) dudit 
signal source avec la somme des vecteurs propres (v^, v'jj 
associes auxdites valeurs propres dominantes, 

- recherche du premier maximum de cette fonction de 
correlation (c(t), f(t)). 

20 2)Procede selon la revendication 1, caracterise en ce que le 
nombre (d, d 1 ) desdites valeurs propres dominantes {X±, X* i) 
est predetermine. 

3) Procede selon la revendication 1 caracterise en ce que le 
rapport de la somme desdites valeurs propres dominantes a la 

25 somme de toutes les valeurs propres est superieur ou egal a 
un nombre predetermine. 

4) Procede selon la revendication 1 caracterise en ce que, 
comprenant de plus une etape d 1 estimation du bruit additif 
(N) dans le canal de transmission, lesdites valeurs propres 

30 dominantes sont telles que leur somme soit inferieure ou 
egale a la somme de toutes les valeurs propres diminuee 
dudit bruit additif (N) . 

5) Procede selon la revendication 4, caracterise en ce que 
1' estimation du bruit additif (N) est realisee en 

35 normalisant le bruit instantane (Nq) qui est lvalue au moyen 
dudit signal de reception ( S ) , dudit signal de reference 
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(TS) et d'une estimation de la reponse impulsionnelle (X) du 
canal de transmission. 

6) Procede selon la revendication 5 caracterise en ce que, en 
notant A la matrice de transmission associee audit signal de 
reference (TS), 1' expression du bruit instantane (Nq) est la 
suivante : N 0 = S - A.X. 

7) Procede selon la revendication 6, caracterise en ce que 
ledit bruit additif (N) est de plus moyenne. 

8) Procede selon l'une quelconque des revendications 1 a 7, 
caracterise en ce que ladite matrice de caracterisation (L) 
resulte d'une operation de lissage. 

9) Procede selon l'une quelconque des revendications 
precedentes caracterise en ce que ledit signal de 
caracterisation est une estimation de la reponse 
impulsionnelle (X) du canal de transmission. 

10) Procede selon l'une quelconque des revendications 1 a 8, 
caracterise en ce que ledit signal de caracterisation est 
ledit signal de reception (S). 
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